Cluster mergers, core oscillations, and cold fronts 
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ABSTRACT 

We use numerical simulations with hydrodynamics to demonstrate that a class of cold fronts in 
galaxy clusters can be attributed to oscillations of the dark matter distribution. The oscillations 
are initiated by the off-axis passage of a low-mass substructure. From the simulations, we derive 
three observable morphological features indicative of oscillations: 1) The existence of compressed 
isophotes; 2) The regions of compression must be alternate (opposite and staggered) and lie on 
an axis passing through the center of the cluster; 3) The gradient of each compression region 
must pass through the center of the cluster. Four of six clusters reported in the literature to 
have cold fronts have morphologies consistent with the presence of oscillations. The clusters with 
oscillations are A496, A1795, A2142, and RX J1720. 1-^2638. Galaxy clusters A2256 and A3667 
are not consistent so the cold fronts are interpreted as group remnants. The oscillations may be 
able to provide sufficient energy to solve the cooling-flow problem and, importantly, provide it 
over an extended duration. 

Subject headings: X-rays: diffuse background — X-rays: galaxies: clusters 



1. Introduction 

The high spatial resolution of the Chandra x- 
ray telescope has led to the discovery of galaxy 
cluster phenomena neither previously seen by x- 
ray telescopes nor previously predicted by numer- 
ical simulations. Among these are cold fronts, 
a term coined in Vikhlinin, Markevitch, & Mur- 
ray (2001) and observed in many clusters includ- 
ing the prototypes Abell 2142 (Markevitch et al. 

2000) , Abell 3667 (Vikhhnin et al. 2001), and RX 
J1720.1-H2638 (Mazzotta et al. 2001) as weh as 
Abell 1795 (Markevitch, Vikhlinin, & Mazzotta 

2001) , Abell 2256 (Sun, Murray, Markevitch, & 
Vikhhnin 2002), and Abell 496 (Dupke & White 
2003). 

Cold fronts share the following two characteris- 
tics: 1) a steep, non-axisymmetric x-ray flux gra- 
dient and 2) a steep temperature gradient with the 
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cooler gas on the denser side in approximate pres- 
sure equilibrium. The density and temperature 
gradients at the front are sufficiently steep to al- 
low them to be modeled as discontinuities, within 
the resolution limits of the detector. 

Cold fronts are not shock fronts. Markevitch 
et al. (2000) estimates M. < V2 for the front in 
A2142. For A3667, Vikhlinin et al. (2001) estimate 
A4 = 1 ± 0.2 from the ratio of the pressure on ei- 
ther side of the front. The same method applied to 
RX J1720.1-h2638 gives M = 0.4l°;^ (Mazzotta 
et al. 2001). 

Steep, non-axisymmetric luminosity gradi- 
ents were observed by the ROSAT x-ray tele- 
scope but misinterpreted as shock fronts (Marke- 
vitch, Sarazin, & Vikhlinin 1999). Indeed, non- 
axisymmetric luminosity gradients had been noted 
in many clusters prior to Chandra (Mohr, Fabri- 
cant, & Geller 1993; Mohr, Evrard, Fabricant, & 
Geller 1995; Buote & Tsai 1996). 

The presently postulated explanation for cold 
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fronts is that they arc the low-cntropy remnant of 
a recently merged substructure (Markevitch et al. 
2000). The low-entropy gas is in the process of, 
or has been stripped via ram pressure from the 
merging substructure. The temperature disconti- 
nuity can arise from either the parent halo suffer- 
ing shock heating, or the low-entropy core cooling 
adiabatically as it expands in response to the loss 
of its potential well. A combination of factors is 
possible. 

The subcluster-remnant scenario begs the ques- 
tion: where is the parent cluster's low-entropy 
core? The cool dense gas of the cold front is the 
remnant of the group. Radiative cooling is more 
efficient in the dense centers of massive clusters. 
This should lead to a large reservoir of low-entropy 
gas. This cool gas is observed in large clusters, 
though not in the abundance anticipated by the 
heuristic argument given which is the classic "cool- 
ing flow" problem. 

This paper proposes cold fronts are a symptom 
of more than one phenomenon. The previously 
postulated explanation that the front is due to gas 
from merging object is valid only for one class of 
fronts which we call here the "remnant class" . We 
suggest another class of fronts in which the front 
is due to oscillations of the core after the merger. 
The significant difference is the cool gas is origi- 
nal, central gas from the parent halo. This class 
is called the "oscillation class" in the rest of this 
paper. This model shares many traits with the 
model proposed in Markevitch et al. (2001). In 
their model, the cool gas is sloshing in the gravi- 
tational potential of the cluster. The difference is 
that in the modc^l prc^sented here, it is the cluster 
gravitational potential that is oscillating, not the 
gas within the potential. The dark matter is mov- 
ing en m,asse in an oscillatory fashion, carrying 
the bulk of the gravitational potential with it. The 
displacement of the gravitational potential induces 
motion in the gas as well. The model of Marke- 
vitch et al. (2001) requires a non-gravitational per- 
turbation to the gas in the core. Such a pertur- 
bation can be provided by an AGN jet or large- 
scale turbulent motion induced by matter in falling 
along a filament (Nagai & Kravtsov 2003). Hence, 
their model potentially represents another class of 
cold fronts. 

Mazzotta et al. (2001) suggests a model which 
is a hybrid of the remnant and oscillation classes: 



a subcluster forms co-spatially with a larger halo 
but offset sufficiently such that the subcluster os- 
cillates within the center of a larger halo. This 
model differs from the remnant class in that the 
cool gas is never stripped from its dark matter halo 
and differs from the oscillation class in that it is 
the subcluster that is oscillating, not the primary. 

There are two advantages to the oscillation 
model presented here. 

1. The oscillations are able to persist long after 
the passage of the merging structure. 

2. Major mergers are not required to instigate 
oscillations. More common high-mass-ratio 
interactions are sufficient. 

Together, these two advantages predict that cold 
fronts should be common, and not unique to clus- 
ters exhibiting recent major merger activity. Ev- 
idence supporting this prediction has been re- 
ported. Markevitch, Vikhlinin, & Forman (2002) 
find 25 of 37 clusters selected from ROSAT vtn&ges 
as relaxed (at least in morphology) have brightness 
edges similar to those in typical cold front clusters. 
Such a large fraction of clusters with cold fronts 
is in conflict with expectations. Fujita, Sarazin, 
Nagashima, & Yano (2002), using an analytical 
model of the mass distribution function, find that 
Va of clusters should have large subhaloes. Not all 
subhaloes should exhibit cold fronts, implying that 
less than Va of clusters should have cold fronts. 

In a set of simulations of preheated galaxy clus- 
ters, Bialek, Evrard, & Mohr (2002) found one 
cluster with a structure reminiscent of cold fronts. 
The cluster was undergoing a recent merger of 
comparably-massed objects, hence belongs to the 
remnant class. They simulated 68 clusters and 
found only one instance of a cold front, suggesting 
that remnant-class cold fronts are rare. 

Nagai & Kravtsov (2003) performed simula- 
tions of cold dark matter cluster mergers specif- 
ically designed to study cold fronts of the rem- 
nant class. They verify the transient nature of 
cold fronts around the substructure remnant. 

Chemical gradients provide a further constraint 
to test the model: in particular, the absence of a 
gradient. The absence of a gradient is difficult to 
explain if the cool, dense gas of the cold front is a 
remnant of an alien structure. However, the oscil- 
lation model predicts no chemical gradient across 
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the cold front unless a chemical gradient existed 
prior to and independently of the perturbing in- 
teraction. The model predicts density gradients in 
regions in which the local gas is compressed. The 
compression increases the gradient of any chemi- 
cal variation that already exists, but does not in 
itself inject gas with a dissimilar chemical history. 
The absence of a gradient in the A496 cold front is 
reported by Dupke & White (2003). Mergers have 
been postulated to destroy chemical gradients (Dc 
Grandi & Molendi 2001); turbulence generated by 
the merger efBciently mixes the gas. However, 
the remnant would necessarily be destroyed by the 
mixing. 

Of the well-studied cold front clusters, we ar- 
gue below that the fronts in A2256 and A3667 
are members of the remnant class, while RX 
J1720.1-F2638, A496, A1795, and A2142 are mem- 
bers of the oscillation class. 

The paper is laid out as follows. In Sec. 2, the 
methodology and initial conditions of the simula- 
tions are described. The results of the simulations 
are described in Sec. 3 and observational conse- 
quences are derived. With the observational crite- 
ria, in Sec. 4 six cold front clusters from the litera- 
ture are examined in detail and classified as either 
remnant or oscillation class. The possibility that 
the energy is sufficient to explain the cooling flow 
problem is explored in Sec. 5. We summarize the 
results in Sec. 6. 

2. Simulations 

We used a series of numerical simulations to ex- 
plore the effects of a merger on the core of a clus- 
ter. Both the collisionless dark matter and the gas 
were evolved using the hydro/N-body code, Hy- 
dra (Couchman, Thomas, & Pearce 1995), which 
models the hydrodynamics using smoothed parti- 
cle hydrodynamics (SPH). The initial clusters and 
groups were extracted from a cosmological simu- 
lation of a (40 Mpc )^ volume. Substructure in the 
haloes was maintained. This provides for a more 
realistic state for the cluster and group haloes. 
The cluster and group were evolved in isolated 
volumes for a few crossing times before combin- 
ing. The short evolution allows artifacts due to 
the abrupt change in the boundary conditions and 
local tidal field to dissipate. 



Simulation A models a high mass-ratio inter- 
action. A group of mass 4.7 x 10^'^ was col- 
lided with a cluster of mass 1.4 x 10^^ for a 
mass ratio of 30:1. The cluster and group were 
extracted from a CDM simulation. The mass res- 
olution of the cluster and group was boosted by a 
factor of two. The cluster was modeled by 371 678 
particles divided 57:43 dark matter to gas, a ratio 
set by evolution within the parent CDM simula- 
tion. The group was modeled by 7940 particles in 
a dark matter to gas ratio of 62:38. The dark mat- 
ter particles had masses 13 times that of the gas 
particles. The group was placed 4.6 Mpc from the 
cluster and given an initial velocity of 500 km s~^ . 
This is a bound state; the bound/unbound tran- 
sition is 1600 km s^^ . The initial impact param- 
eter (defined as the separation vector of closest 
approach in the absence of gravitational accelera- 
tion) was 1 Mpc, leading to a closest approach of 
110 kpc. There was no radiative cooling in Sim- 
ulation A. The cluster-group system was evolved 
for an elapsed simulation time of 4 Gyr . 

To boost the resolution of Simulation A, each 
particle was split into two and distributed ran- 
domly within a local radius in which there were 
originally 32 particles. For dark matter particles, 
the particles were given random velocities such 
that the total momentum, kinetic energy, and lo- 
cal velocity dispersion (determined from the 32 
closest particles) of the original particle were con- 
served. 

Simulation B models an intermediate mass- 
ratio interaction. A group of mass 4.3 x 10^"^ Mq 
was collided with a cluster of mass 4.2 x 10^'^ M© 
for a mass ration of 10:1. The cluster and group 
were extracted from a A-CDM simulation (O^ = 
0.3, A = 0.7). The simulation used 233 461 parti- 
cles split almost 50/50 into dark matter and gas. 
The ratio of the masses of the dark matter to gas 
was 5.7:1. The group was placed 6.2 Mpc from the 
cluster and given an initial velocity of 500 km s~^ . 
This is a bound state. The bound/unbound tran- 
sition is 763 km s^^ . The initial impact parameter 
was 1.5 Mpc, leading to a closest approach of 356 
kpc. In contrast to Simulation A, radiative cool- 
ing was included. The cluster-group system was 
evolved for an elapsed simulation time of 20 Gyr . 
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3. Simulation results 

Oscillations in the core are induced by the pas- 
sage of the group. The core moves with respect to 
the outer halo which, in turn, begins to move with 
respect to the core. The details arc illustrated 
in Fig. 1 which plots contours of equal projected 
mass at time intervals separated by 50 Myr from 
Simulation A. The group approaches from the left 
in frames 1 and 2 (group not seen), passes over the 
cluster and is deflected to the lower right where it 
continues out of the frames. Closest approach is 
110 kpc in frame 3. Compression of the contours 
are seen in frames 4, 6, and 7. In frame 4, the 
compression is seen on the upper face, 16 kpc in 
the direction of the impact radius and is due to 
the movement of the cluster center toward the ini- 
tial impact point. In frame 6, the cluster center is 
now stationary, but the halo has moved upward, 
overshooting the core's displacement and causing 
compression on the lower side 16 kpc away. The 
overshoot re- accelerates the core upward, causing 
the core to compress material above itself, 8 kpc 
away, while the halo is still moving upward, con- 
tinuing the compression in the lower half at the 
larger radius of 32 kpc. Compare frame 7 with 
Fig. 7 where there are not two but three regions 
of compression laid out in a similar fashion. By 
frame 8, 300 Myr after the initial impact in frame 
2, the core and halo have come to rest centered on 
a point displaced 50 kpc above their initial posi- 
tion. Note that the core oscillates perpendicularly 
to the initial velocity vector of the impactor. 

At this point it is important to qualify what 
we are calling "oscillation". No single iso density 
contour oscillates in a periodic fashion (like a pen- 
dulum) with respect to the fixed reference frame. 
Instead, motions are periodic only with respect to 
the cluster center, which is not a fixed reference 
frame. 

From Fig. 1, we estimate the duration of the 
oscillations is 250 Myr for Simulation A. For the 

smaller mass ratio in Simulation B, the oscillations 
persist for at least 600 Myr but have dissipated by 
1.2 Gyr. 

Crucial to the model presented here are the fol- 
lowing three observational consequences: 

1. There can be one, two, or three regions of 
compressed isodensity contours. 



2. The regions of compression must be alter- 
nate (opposite and staggered) of each other 
and lie on an axis which passes through the 
center of the cluster. 

3. The gradients of compressed isodensity must 
point through the center of the cluster. 

The remnant model cannot explain more than one 
region of compressed isodensity contours. Coinci- 
dental mergers or separate substructures can pro- 
duce two cold fronts. However, they need not be 
aligned along any preferred axis or be alternate. 
The model presented by Markevitch et al. (2001), 
in which the gas is oscillating in a fixed potential 
set by the dark matter, predicts the third conse- 
quence, but does not explain the ringing which 
leads to the first two consequences. Indeed, it is 
not obvious that a model in which the dark mat- 
ter at the center is oscillating within a global po- 
tential, dragging the gas with it, will reproduce 
the ringing. The alternate compression regions are 
due to the delayed response of the outer halo to the 
perturbation of the halo by the passing substruc- 
ture. We will return to these three consequences 
in Sec. 4. 

There exists the opportunity for confusion if 
the cluster-substructure interaction induces oscil- 
lation of the core and the substructure remains in, 
or subsequently falls back on the core to produce 
both a sloshing and merger-remnant cold front. 
Such a cluster would have an merger-remnant cold 
front while the core of the primary cluster would 
be undergoing oscillations. But these two cold- 
front phenomena are still distinct, even if they can 
exist in the same cluster at the same time. An 
oscillation in a cluster with the merger remnant 
in the immediate vicinity has not been explored 
by these simulations. Given that the presence 
and strength of the oscillation-type cold fronts is 
a function of time, it is likely that the dominant 
cold front class in a mixed scenario would be very 
sensitive to the stage of the merger. 

There are other observational consequences 
which are not discussed here: isophotal twisting 
and multiple oscillation modes. The possibility 
of isophotal twisting occurring exists, but is not 
explored in this paper since the effect is not as 
obvious. However, we suspect that it would be a 
further probe of the merger geometry. Multiple 
sets of oscillations would be initiated by multiple 
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merger events (either separate objects or a single 
object on separate passes). If multiple oscillation 
modes are concurrent, the resultant morphology 
will be distorted making analysis problematic. 

The core velocity relative to the simulation vol- 
ume for Simulation B is illustrated in Fig. 2. The 
strongest perturbation to the core of the primary 
halo is in the direction of the impact parameter, 
perpendicular to the initial group velocity vector. 
The maximum velocity is 300 km s^^ which is less 
than the speed of sound of the cluster medium 
(c^ ~ 340km s-i). 

Fig. 3 illustrates the displacement of isoden- 
sity contours in response to the perturbation, with 
respect to the center of the primary cluster halo 
for Simulation B. Plotted are the fractional shifts 
[dr/r) in the centers (along the impact radius vec- 
tor) of isodensity surfaces of varying radii at dif- 
ferent times. Closest approach of the perturbing 
group is 300 kpc. Non-axisymmetric density gra- 
dients are created where the slopes of the curves 
are steepest. Negative slopes indicate isodensity 
compression, positive slopes indicate isodensity 
stretching. Before impact (600 Myr prior to clos- 
est approach), the maximum fractional variation 
is < 10%. At closest approach, there is one large 
displacement of ~ 60% for the isodensity contours 
at a radius of 20 kpc. At 600 Myr after closest ap- 
proach, there are three regions of isodensity con- 
tour compression: at 8 kpc, 35 kpc, and at > 80 
kpc. By 1.2 Gyr, there is only a compressed gra- 
dient left in the center (6 to 8 kpc) and at large 
radii (> 200 kpc). 

The relative position of the gas and the dark 
matter haloes clarifies if the gas is simply being 
dragged along with the dark matter or is also be- 
ing displaced. Figure 4 maps the ratio of the dark 
matter mass to gas mass at a time 200 Myr af- 
ter closest approach in Simulation A (frame 7 of 
Fig. 1). The enhancement of the dark:gas ratio 
is evidence that it is the dark matter distribution 
that is oscillating, bringing the gas with it. How- 
ever, there is a lag between the displacement of the 
dark matter and the displacement of the gas. The 
compression of the gas on the forward side leads to 
a gas pressure force which displaces the gas from 
the center of the potential well. The displacement 
leads to adiabatic cooling of the gas on the lee- 
ward side and subsequent as the gas is pushed out 
of the potential well. 
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Fig. 1. — Oscillations induced by the passage of a 
merging sub-cluster in Simulation A. The contour 
lines map projected gas mass. Contour lines are 
spaced by a factor of V2. The dotted lines give 
a constant frame of reference. Frames are 50 Myr 
apart. The scale of each box is 200 kpc. 



350 p 
300 
250 
200 
150- 

cn 

E 100 
50 

o| 

-50 

-100- 
-150- 



X (initial clump velocity) 
y (impact parameter vector) 
z (normal) 




10 

time (Gyr) 



15 



20 
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est approach). Only the shift in the direction of 
the impact parameter is measured. For Simulation 
B. 




Fig. 4. — Mass ratio map in Simulation A. This 
frame corresponds to frame 7 in Fig. 1. 



The sharp discontinuity of x-ray emission ob- 
served in cold-front clusters is an important fea- 
ture of cold fronts that any theory must predict 
and explain. The discontinuity is entirely at- 
tributable to a change in the gas density across 
the front, with the higher-density gas on the side 
of the front with the higher emission (contrary to 
the situation in which the emission excess is due 
to the presence of a shock) . 

The simulations have yet to recreate the sharp- 
ness observed. Figure 5 illustrates the density pro- 
files across the compression regions in a simulated 
cold front. Two profiles are plotted, derived from 
the gas particles lying in conical regions extend- 
ing outward and opposite from the center of the 
cluster. The cluster is as in frame 6 in Fig. 1 with 
the cones having opening widths of 20°. The two 
cones extend above and below the cluster to con- 
trast the density profiles in these regions. 

At the resolution of the simulations, the simu- 
lated fronts are consistent with the sharpness of 
the discontinuities observed. The simulation spa- 
tial resolution is provided by the SPH method. 
The SPH smoothing length is set by the local num- 
ber density of particles. The mean SPH smoothing 
radii for the particles as a function of radius is il- 
lustrated by the lines at the top of Fig. 5. The 
mean SPH smoothing lengths in these regions are 
comparable to the width of the compression re- 
gions in the simulations. 

It is as yet unknown whether increased spa- 
tial resolution provided by simulations with larger 
number of particles will lead to sharper density 
discontinuities. According to the model presented 
in this paper, higher resolution simulations should 
produce sharper fronts. The absence of sharper 
fronts in higher-resolution simulations will falsify 
the model. 

4. Comparison with observations 

The observational consequences of the model 
we have proposed permits morphological-based 
classification of cold-front clusters into oscillation 
and non-oscillation classes. More than morpholog- 
ical data can be used in such a classification. The 
velocity of the cold front is, on its own, a better 
discriminant. However, determining the speed of 
a cold front from the x-ray data is not a trivial 
exercise. It requires high quality data with spec- 
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troscopic information and assumptions about the 
relation between the temperature-density state of 
the gas and the velocity of the front. The simple 
morphological rules described in this paper may 
be applied directly to the image data. 

In the literature, six clusters are described 
in sufficient detail to permit comparison with 
the model proposed here: A496, A1795, A2142, 
A2256, A3667, and RX J1720.1-I-2638. We ex- 
amine each of them here. Fig. 6 shows the x-ray 
flux for the six clusters. The data are from, re- 
spectively, ACIS-S3 chip over two observations, 
ACIS-S3 chip over three observations, ACIS-S3 
chip over two observations, ACIS-I array (1 obser- 
vation) and ACIS-S3 chip (3 observations), ACIS- 
I array over two observations, and ACIS-I array 
over one observation. The images were created 
with an adaptive smoothing kernel. For 5 of the 
clusters, Reiprich & Bohringer (2002) have pub- 
lished masses: A496 (5.1 x lO^ftes^Me), A1795 
(11.9X 10"/i-i M0 ), A2142 (14.7 x lO^fteg^ Mq ), 
A2256 14.9 x IQi^/ies^ Mq ), and A3667 (9.6 x 
IO^^Hq^ Mq ). The range in cluster masses is sim- 
ilar to the range in the simulations. 

A496 (Dupke & White 2003) (Fig. 6) has a 
pair of non-axisymmetric isophotes discontimiitics 
which are alternate and lie along an axis through 
the center of the cluster. The isophotes are com- 
pressed Uh^rl kpc (1") SSE and 66h^^^ kpc (1.5") 
NNW of the cluster center. All three observational 
consequences of the oscillation model are satisfied. 
The model for the remnant class of objects, in 
which the group remnant is passing through the 
primary clusters halo, predicts a mctallicity dis- 
continuity at the front. But the model for the 
oscillation class is inconsistent with a mctallicity 
discontinuity since, unless there already exists an 
axisymmetric gradient, it is the in situ gas that 
is being compressed. No mctallicity discontinu- 
ity was found across the cold front in A496 by 
Dupke & White (2003). We classify A496 as an 
oscillation- type cold front. 

A1795 (Markevitch et al. 2001) (Fig. 6) has a 
broad front at its southern edge and a tighter front 
(Fig. 7) at the northern side of the core in what 
is described as a filament. The isophotes are com- 
pressed 11/i^^kpc (0.13") N and Sih^^kpc (1") 
S of the cluster center. The fronts lie along the 
north-south axis. Hence, the three observational 
requirements for an oscillation class cold front are 



satisfied and we classify A1795 as such. 

A2142 (Markevitch et al. 2000), x-ray fltix map 
shown in Fig. 6, exhibits a clear compression of its 
isophotes Slh^l kpc (0.7") SE of the cluster cen- 
ter and at larger radii SOO/igg^kpc (2.5") NW of 
the center. The compression regions are alternate 
and lie along an axis passing through the clus- 
ter center. We classify A2142 as containing an 
oscillation-class cold front. 

A2256 (Sun et al. 2002) (Fig. 6) is revealed in 
the x-ray band to contain an irregularly shaped 
knot of cool dense gas located to the west of the 
cluster center. The sharpest edge of the cool gas 
faces SSE. There is only one edge. The gradient of 
the edge is not along an axis through the cluster 
center. A2256 does not contain an oscillation-type 
cold front. We classify A2256 as a merger-class 
front. 

A3667 (Markevitch et al. 1999; Vikhlinin et al. 

2001; Mazzotta, Fusco-Fcmiano, & Vikhlinin 
2002), shown in Fig. 6, has a single cold front 
410/ig5^ kpc (5.5") SE of the cluster center. It also 
contains a filament of gas to the NE of the cluster 
center. Though the gradient in the region of the 
front does point through the cluster center, the 
single front and the filament are more consistent 
with an ongoing major merger. Simulations re- 
produce the morphology with a massive head-on 
merger (Roettiger, Burns, & Stone 1999). We 
classify A3667 as a merger-class front. 

RX J1720.1-F2638 (Mazzotta et al. 2001) is an 
almost relaxed cluster (see Fig. 6). There is a clear 
compression of the isophotes at three locations: 
two to the SE at Uh^^ kpc (0.2') and 240h^^ kpc 
(1.1') from the center and one SShg^kpc (0.4') to 
the NW. All are along a position angle of 140°, 
compressed perpendicular to the position angle. 
The position and alignment of the compression re- 
gions is consistent with oscillation-class fronts. 

Of the six clusters reported in the literature to 
have cold fronts, four (A496, RX J1720. 1+2638, 
A2142, and A1795) are consistent in morphology 
with oscillation-class cold fronts. Two are incon- 
sistent (A2256 and A3667) and we classify with- 
OTit further examination as having remnant-class 
cold fronts. Six clusters are too few to permit 
a satisfying statistical sample. However, if ^3 of 
clusters have cold fronts (Markevitch et al. 2002), 
while less than ^/s are expected to have merger 
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remnants within their haloes (Fujita et aL 2002), 
then more than half the cold front clusters must 
originate from some process other than remnant 
group haloes. 

In the simulations, the compression regions 
were found over a range of scales from 8 to 80 kpc. 
The observed compression regions in the clusters 
classified as oscillation type were over a similar 
range in magnitude but overall larger scale: 11 
to 300/ig5^ kpc . The smallest oscillation feature in 
the simulations (8 kpc) would correspond to only 
10" in the Chandra observation of the nearest clus- 
ter (A496). Chandra is not limited by resolution 
in resolving such small features, but it is limited 
by counts. Fig. 8 reveals the inner 3' region of 
A496. There is a region of compressed isophotes 
15/ig5^ kpc (21") E of the x-ray peak. The region 
is not collinear with the other two regions of com- 
pression to the north and south of the cluster cen- 
ter, but the gradient vector is parallel with a line 
passing through the x-ray peak, raising the pos- 
sibility that A496 is undergoing two concurrent 
oscillations. 

5. Energy in the oscillations and cooling 
flows 

The temperature of the gas in the centers of 
cooling flow clusters are consistently measured to 
not fall below ^/sTyiriai (Fabian, Allen, Crawford, 
Johnstone, Morris, Sanders, & Schmidt 2002). 
This runs contrary to expectation; the radiative 
efficiency of cool dense gas does not fall off un- 
til below IC^K. To prevent gas from cool- 
ing catastrophically requires an energy injection 
method that is both plentiful and not instanta- 
neous. Fabian et al. (2002) notes a required rate 
of energy injection of up to 10^^ erg sec~^. The 
kinetic energy in the oscillations can both provide 
a source of energy and spread this energy injection 
over a period of time. But is it enough? 

An examination of the energy stored as kinetic 
energy in the three phases demonstrates that os- 
cillations could be the unknown source of power 
that regulates the central temperature in cooling 
flow clusters provided there exists a mechanism 
to transfer kinetic energy from the galaxies. The 
cluster in Simulation B has a mass of 4.2 x 10^^ Mq 
which is not particularly massive. Clusters with 
similar masses are observed to have luminosities 



of about 10'*'*erg sec~^. Fig. 9 graphs the kinetic 
energy stored in the dark matter, galaxies, and 
gas in Simulation B. Closest approach occurs at 
7.2 Gyr. For the gas, a minimum of 4 x 10'"'^ erg 
is stored as kinetic energy in the gas phase. The 
gas disperses the energy in 700 Myr for a mean 
heating rate of 1.8 x 10'*'^ erg sec~^. This is one 
fifth of the required injection rate. However, the 
other two phases: galaxies and dark matter, store 
substantially more energy. The relative motions 
of galaxies store a minimum kinetic energy of 
1 x 10^^ erg and disperse this energy over a longer 
time-scale, ^ 2.4 Gyr, for a mean injection rate 
of 1.3 X lO**** erg sec~^ which is of the same order 
of magnitude as the luminosity. The largest reser- 
voir of kinetic energy is the dark matter. The dark 
matter receives at least 3.6 x 10^^ erg from the in- 
teraction which it dissipates in 1.2 Gyr for a mean 
dissipation rate of 1.0 x 10'*''' erg sec""'^. Though 
this is the well above the amount of power re- 
quired, there are no efficient means for transferring 
energy between the two phases. If a mechanism 
can be found to efficiently transfer galactic kinetic 
energy into gas thermal energy, the energetics of 
the system are sufficient to maintain the minimiim 
temperature in the centers of cooling flow clusters. 

6. Summary 

The cluster morphology observed in some clus- 
ters with cold fronts has been shown to arise from 
oscillations of the cluster core induced by gravi- 
tational interactions with subclusters. Numerical 
simulations with hydrodynamics give us a list of 
morphological signatures of these oscillations: 

1. There can be one, two, or three regions of 
compressed isodensity contours. 

2. The regions of compression must be alter- 
nate of each other and lie on an axis which 
passes through the center of the cluster. 

3. The gradients of compressed isodensity must 
point through the center of the cluster. 

We associate the compressed isodensities with cold 
fronts and show that 4 of 6 clusters with cold fronts 
examined exhibit these signatures. The other two 
clusters which do not exhibit these features are 
likely undergoing a merger of a subcluster with 
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the cold gas coming from the remnant of the sub- 
cluster. 

That the oscillations create cold-front mor- 
phologies solves the discrepancy between the num- 
ber of clusters with cold-fronts observed and the 
number expected if cold-fronts were due only to 
the remnants of merged subclusters. The oscilla- 
tion class has two major advantages to explaining 
the increased frequency: 

1. The oscillations are able to persist long after 
the passage of the merging structure. 

2. Major mergers are not required to instigate 
oscillations. More common high-mass-ratio 
interactions are sufficient. 

The energy contained in the kinetic energy of 
the oscillating intracluster medium is insufficient 
to offset that lost by radiative cooling and hence 
does not solve the cooling flow problem. How- 
ever, the kinetic energy carried by the galaxies in 
the oscillation is of the correct magnitude. In ei- 
ther case, the duration over which the energy is re- 
leased to the cluster intracluster medium is ~ 100 
Myr to 1000 Myr, providing a steady source of 
energy. 
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Fig. 5. — Density profiles for gas on opposite sides 
of the cluster center. See the text for descriptions 
of the cones. The mean gas smoothing radii are 
denoted by the horizontal lines. The first row is 
for the top cone. The second row for the bottom 
cone. This frame corresponds to frame 6 in Fig. 1. 
The density profiles do not have the sharp discon- 
tinuities the observations suggest, but the profiles 
are consistent with abrupt discontinuities softened 
by the SPH smoothing length. 




Fig. 6.— Galaxy clusters A496, A1795, A2142, 
A2256, A3667, and RXJ1720. 1+2638. The scales 
are, respectively, 1' = 44/ig5^ kpc , 1' = SAh^^ kpc , 
1' = 122/1^:/ kpc, 1' = 78/i65^kpc, 1' = 
75/ig5^ kpc , and 1' = 220/ig5^ kpc . 



Galaxy Cluster A1 795: Chandra ACIS: Flux, 0.6 - 7 keV ^ „-3 




13:49:00 13:48:55 13:48:50 

Fig. 7. — A detail of galaxy cluster A1795. 
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Galaxy Cluster A496: Chandra ACIS: Flux, 0.6 - 7 keV ^ „-3 




4:33:45 4:33:40 4:33:35 

Fig. 8. — Detail of galaxy cluster A496. 




Fig. 9. — Kinetic energy in the dark matter, galax- 
ies, and gas for Simulation B. 



